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I. INTRODUCTION 

Charge-transfer reactions involving few particles in atomic physics are very challenging 
and interesting from both theoretical and experimental points of view and we study here the 
problem of charge transfer in some atomic reactions involving exotic particles. Specifically, 
we study muon transfer in D— H^, T— H^, and T— D^ systems where the suffix /i denotes 
a muonic atom with the electron replaced by a muon (/i~). We also study the problem of 
antihydrohen (H) formation in antiproton-positronium collision with the positronium (Ps) 
in an initial Is state. 

On the theoretical side, in these transfer reactions one needs to consider rearrangement 
of a charged particle. Because of the Coulomb interaction one needs a careful treatment 
of the dynamics for a correct description. If one can identify the basic dynamical ingredi- 
ents necessary for a satisfactory description of these processes involving a small number of 
particles, such a study will help us to formulate models in more complex situations. On 
the experimental side, the present study involving muon and positron transfer is of current 
interest in the muon-catalyzed fusion cycle [|l]-§] and in the formation and study of the 
antihydrogen atom [^§. 

Although there are some experimental measurements and several theoretical investiga- 
tions on these processes, there still remain discrepancies among various studies. Here we use 
a different theoretical approach based on a detailed few-body dynamical consideration for a 
careful reinvestigation of these three-body charge-transfer reactions. Traditionaly, such prob- 
lems are investigated by a tractable approximation scheme in the Schrodinger framework, 
without explicitly considering a few-body dynamical equation. In addition to variational 
calculations, these schemes include close-coupling, hyperspherical, and adiabatic approxi- 
mations. Here we would like to point out that the processes of muonic transfer reactions 
and antihydrogen formation are three-body Coulombic rearrangement collisions. Conse- 
quently, it seems reasonable that in addition to approximations based on the Schrodinger 
equation, a detailed few-body consideration is useful. In what follows we develop a method, 



which is based on detailed few-body equations rather than the effective potential treatment 
employed in alternative investigations. 

For the three-charged particle system, say TDyU, only two asymptotic configurations are 
possible, i.e. (Dyu) T and (Tyu) D. This suggests to write down a set of two coupled equations 



for components \l/i and ^2 of the wave function \1' = \E'i + \E'2 PJIHj] with each component 
carrying the asymptotic boundary condition for a specific configuration. One such equation 
with two components for the three-particle system was first written by Hahn |^ following 
the most general decomposition of the three-body wave function into three components 
suggested by Faddeev ||ll[ and is usually referred to as the Faddeev-Hahn equation [|l^]. We 



solve the integro-differential form of this equation by a six-state close-coupling approximation 
scheme which consists in expanding the wave function components \I'i and \l/2 in terms of 
eigenfunctions of subsystem Hamiltonians in initial and final channels, respectively. The 
resultant coupled equation is then projected on the expansion functions. After a partial- 
wave projection this leads to a set of one- dimensional coupled equations for the expansion 
coefficients, which is solved numerically. 

Recently, there have been considerable theoretical and experimental interests in the study 
of the muon-transfer reactions between hydrogen isotopes in the muon catalyzed fusion cycle 

(1) 
^ T^ + D . 

The measurements for the transfer rates 

^tr = ^tr^^o, (2) 

with (T^f being the transfer cross section, v the relative velocity of the incident particles and 
No = 4.25 X lO^^cm"^ the liquid hydrogen density, are listed in Table I together with recent 
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T + D^- 


-^T^ + D. 



theoretical calculations. One can see differences between different experimental data [p!^ -[T9| 



and theoretical results 20-25 . 



One of the most attractive reactions for H formation is the three-body positron-transfer 
process 

p + Ps^H + e". (3) 

Although no experimental cross sections are available, this process is being used at CERN 
for the production and study of antihydrogen. A number of calculations have recently been 
carried out to calculate the cross section of reaction (§) as a function of the incident Ps 
energy. The calculations were performed by different methods, for instance, with hyper- 
spherical coupled-channel expansions p6| and close coupling approximations (CCA) 



As an additional test of the present method, calculations for the S— wave antihydrogen 
formation (^) at low energies are also performed. 

In Sec. II we develop the formalism. The results obtained for reactions (|l]) and (^ are 
given in Sec. III. Finally, we present some concluding remarks in Sec. IV. 

II. THEORETICAL FORMULATION 

Let us take the system of units to be e = ^ = m^ = 1 and denote, say T by 1, D by 2 
and muon by 3. Below the three-body breakup threshold, following two-cluster asymptotic 
configurations are possible in the system (123): (23) — 1 and (13) — 2. These configurations, 
denoted simply by 1 and 2, respectively, are determined by the Jacobi coordinates ifjz.pk) 

rj3 = rg - rj, pk = (rg + mjfj)/{l + m^) - f/,, j ^ k = 1, 2, (4) 

Tj, rrij are coordinates and masses of the particles j = 1, 2, 3, respectively. 

Let us introduce the total three-body wave function as a sum of two components 

^(ri,r2,r3) = ^i(f23,pi) + ^2(^13,^2), (5) 

where "^1(^^23, Pi) is quadratically integrable over the variable r23, and "^2(^^13, P2) over the 
variable ri3. To define \1// (/ = 1, 2) a set of two coupled equations can be written down 



{E~Ho- 1/23)^1(^23, Pi) = (V"23 + 142)^2(^13, P2) 

(6) 

{E-Ho- ri3)^2(ri3,P2) = (^13 + V'l2)^l(r23,Pl) , 

where E is the center-of-mass energy, Hq is the total kinetic energy operator, and Vij{rij) 
are pair-interaction potentials {i y^ j = 1,2,3). Equations (|^) satisfy the Schrodinger 
equation exactly and for energies below the three-body breakup threshold they possess the 
same advantages as the Faddeev equations, since they are formulated for the wave function 
components with correct physical asymptotes. 

In the general case a component of the three-body wave function has the asymptotic form 
which includes all open channels: elastic/inelastic, transfer and breakup. In this case each 
component of the total wave function carries a specific asymtotic behavior. The component 
\I'i carries the asymptotic behavior in elastic and inelastic channels: 

^i(r23,Pi) - e<>i(f23) + E^n^''(^Pi)e^'"'''VPl^n(r23). (7) 

The component "^2 carries the asymptotic behavior in the transfer channels: 

^2(ri3,p'2) -_E^:^(^P2)e^'-^VP2V^™(ri3), (8) 

m 

where e^^^ ^V'i(^23) is the incident wave, ipnlfj^) the n-th excited bound-state wave function 



of pair (j3), A;W = ^2Mi{E - E^J^) with M"^ = m'^ + {l + mj)'^. Here E^-^^ is the binding 
energy of (j3), i 7^ j = 1,2, y4'^'/™(fipj and A^^iVLp,^) are the scattering amplitudes in the 
elastic/inelastic and transfer channels. This approach simplifies the solution procedure and 
simultaneously provide a correct asymptotic behaviour of the solution below the 3-body 
breakup threshold. 

Let us write down Eqs. (|^) in terms of the adopted notations 

^ + ^ + "^ " ^^'^ "^'^^^'^ P'^ = ^^^'^ + ^.^)^.(^l3, Pi) , (9) 

here j 7^ A; = 1, 2 and Mj^^ = n\^ + (1 + f^j)~^ , P/^ = 1 + i^J'^- 



For solving Eq. (^ we expand the wave function components in terms of bound states in 
initial and final channels, and project this equation on these bound states. This prescription 
is similar to that adopted in the close-coupling approximation. Specifically, we use the 
following partial-wave expansion 



^fc(r,3,Pfc)= E ^fMXi{pk.r,:i){Y^{h)®yi{rjz)}LM ' (10) 



LMXl 



m'm 

where C's are the Clebsch- Gordon coefficients, and Y^s are the usual spherical harmonics, 
and L,X,l and M,m',m are the appropriate angular momenta variables and their projec- 
tions. Next we make the following close-coupling-type approximation for the radial part in 
terms of the bound-state wave functions in the initial and final channels: 

^iLiiPk, r,3) - - E fiT\pk)R^:hr,3) , (12) 

Pk n 

where radial components of the bound-state wave functions W^i (rjs) satisfy 



'Ei'^^ + ^ 



^ -ir%^) - 1(1 + I) 



y{k) 



Vj3\Rl:i'{rj3) = 0. (13) 



^Pjr% [dr,s' ''dr^s' 

Then we substitute Eqs. (|TU|)-(|r^ into Eq. (|^), multiply the resultant equation by 
the appropriate biharmonic functions and the corresponding radial functions R^i (vj^), and 
integrate over the corresponding angular coordinates of the vectors fj^ and pk- Then we 
obtain a set of integral differential equations for the unknown functions fni\{pk) 



2M,(E - E^^)f^\p,) + (1^ - ^^^]fj:\pk) = 2M,Y.rdr, 



wh 



drrs dpARl^i\r,3) {Fa(p.) ® V'K^.s)}^, (^.3 + V,k) {Yyip,) ® l/'(ru)}iM 
J pj 

Rl!'\'irk3)fihPj) ■ (14) 

For brevity we have defined a = nlX and a' = n'l'X', and omit the conserved total angular 
momentum label LM. The functions fa\pk) depend on the scalar argument, but Eq. (|T3p 
is not yet one-dimensional. We are using the Jacobi coordinates 



Pj = rj3 - Pkrk3, rj3 = -{PkPk + Pj), rjk = -{ajpj - akPk) , (15) 

with 

Pk = j^^, ak = l-(3k, 7 = l-A/3i, Jy^k = l,2. (16) 

l + ruk 

This shows that modulus of pj depend on two vectors: pj = •yfj^ — PkPk- The integration in 
the right-hand side of Eq. (0) is done over these two vectors. 

To obtain one- dimensional integral differential equations, corresponding to Eq. (0), we 
proceed with the integration over variables {pj,Pk}, rather than {rj3,Pk}- The Jacobian of 
this transformation is 7^^. Thus, we come to a set of one-dimensional integral differential 
equations 

2M,iE - E^^)f^:\p,) + (i; - ^i^l /f (p,) = 

I '^Pk Pk ) 

^T.rdp,S^^{Pk,P^f^{p^. (17) 

where functions S^^^, {pk, Pj) are defined as follows 

Si'^){Pk, P,) = 2pkP, J dp^ J dp,Ri]\r,s) {Yxipk) Yi{f,,)}l^ (^,3 + V,k) 

X {Yyip,) ® Ylifk3)},^M Rn'l'irk3) ■ (18) 

The fourfold multiple integration in equations ( ]T8| ) leads to a singlefold integral and the 
expression (|T8|) for any value orbital momentum L becomes 

sL'^)iPk,P,) = ^^[(2A + 1)(2A' + l)]'^pkp, £dusmujRi\\r,3)iV,3ir,3) + V,kir,k)) 

Rli^i'irk3) E ^w(0, ^, 0)C,^oT^Cf,^;^,y,^(z/„ 7T)Y;^,iuk, n) , (19) 

mm' 

where D^^,{0, u, 0) are Wigner functions, u is angle between pj and pk, Vj between r^s and 
Pj, z/fc between r^s and pk- 

Finally, the set of integro-differential equations for the unknown functions fniliPk) can 
be written as 






dp: 



f^iPi) = g^i: 



dpi'fa'\Pi') / du sin uRl^}{\fi, si) 
"'0 [ |ri/3| \rii 

p^Pi' E ^w(o,^,o)c,^oTm^f;-;„,y,:;,(z/„7r)yzw(^.',vr) 



1 1 

1^ + ^- 



(2A+1)(2A' + 1) 
2L + 1 



(20) 



Here i 7^ i' = 1, 2, gf^ = 47rMi/73, A;W = ^j2Mi{E - eP), uj is angle between the Jacobi 
coordinates pi and /?«', i/j is the angle between rj/3 and pi, i/j' is angle between rj3 and pi' 
with 



sm Vj 



smcj 



and 



cos z/,- 



-(APi + Pi'COSu;). 



(21) 



To find unique solution to system (pOD, appropriate boundary conditions are to be con- 
sidered. First we impose /„; (0) = 0. For the present scattering problem with 1 + (23) 
as the initial state, in the asymptotic region two solutions to Eq.(|20|) satisfy the following 
boundary conditions 



<•(!) 



fls\Pl 



c{2)i 



.(1) 



(1) 



pi-^+00 



sin(A;i pi) + Ku cos{kl pi 



,(2). 



(22) 



fls {P2) ~ Jvi/v2Ki2COs{k\'p2) 



where 1 refer to channel 1 + (23), 2 to channel 2 + (13) and K denotes the corresponding 



on-shell K- matrix ||28|. For scattering with 2 + (13) as the initial state, we have the following 
conditions 

f fil\p.) 



Ji) 



~ JV2/V1K21 cos{k\ 'pi] 

pi— ++CO ^ 



fifip. 



,{2) 



sm{k\ ' P2) + K22 cos{k[ V2 



(2), 



(23) 



P2^+00 

where Vi, i = 1,2 are velocities in channel i. With the following change of variables in 
Eqs.O) 



f(i) 



Ki), 



M 



(2), 



flr(pi) = /iV(Pi) - Mk\"pi) and frj{p2) = nt>{p2) - sin{kY>p2 



c(2). 



,(2), 



(24) 



we can obtain two sets of inhomogeneous equations which are solved numerically. The cross 
sections are given by 
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An 
a. 



'^ k^'^^ 



K 



1-iK 



i^ 'mP1±3 (25) 



k\'>HD-ir + {Kn + K^r-' 



where i = j = 1,2 refer to the two channels and 

D = det K = Kn K22 - K12K21 . (26) 

When kl ^0: K12 = K21 ^ ki\ Kn ~ k\\ in this case a^j. = (T12 ~ l/Zci , and 
Cgi = o"ii ~ const. For comparison with experimental low-energy data it is very useful to 
calculate the transfer rates @) because A^j^(A;| ^ ^ 0) ~ const. 

III. NUMERICAL RESULTS 

To solve the integro-differential equation, one has to calculate the angle integrals in Eq. 



(pop which are independent of the energy E. One needs to calculate them only once and 
store on hard disk for the calculation of other observables, for instance, the cross sections 
at different energies. Subintegrals in Eq. (pO| ) have strong dependence on pi and p^/ {i 7^ 
i' = 1, 2). To calculate S^lJ{pi, pii) at different coordinates an adaptable algorithm has been 
used. In this case using the relation 

x'^ - Pi pi - pI 



cosuo 



(27) 



2APiPi' 

the angle dependent part of Eq. (20) can be written as the following integral 
(„.) _ 4vr [(2A + l)(2A- + l)]i ^^p.+p. (, [ _j^ 

R'niinzix)) E ^w(0, c^(x), Q)C^x^^Cl:J^vm'yL{^r{x), 7:)Y,^,{v,{x) , vr). (28) 

mm' 

Note that the expression ( p8|) differs from zero only in a narrow strip when pi ^ p,/. 

We employ muonic atomic unit: distances are measured in units of a^, where a^ is 
the radius of muonic hydrogen atom. The integro-differential equations were solved by 
usual numerical procedure by discretizing them into a linear system of equations, which are 
subsequently solved by Gauss elimination method. In solving these equations distances upto 



9 



50a^ were considered and 400 — 600 points were used in the discretization. The following 
mass values are used in the unit of electron mass: uih = 1836.152, mo = 3670.481, uit = 
5496.918 and the muon mass is rra^ = 206.769. 

Tables II, III, and IV include our results for the muonic transfer cross sections and rates 
for all hydrogen isotopes (^ using different approximation schemes. We present results for 
two-, four-, and six-state approximation where we include Is, ls+2s and ls+2s+2p states 
of the muonic atoms in the initial and final channels, respectively. In solving the equations 
we employed only the lowest partial wave, e.g., L = 0. As we shall mainly be concerned 
with the experimental muon transfer rates at very low energies, the higher partial waves are 
expected to have neghgible contribution. The 2p states are found to contribute significantly 
in T-D^, moderately in D-H^, and little in T-H^ systems. This is in agreement with similar 
conclusion of Ref. [Q in the T-D^ system. This could be understood qualitatively from the 



following consideration. At zero incident energy the relative velocity in the final state after 
muon transfer is the highest in the case of T-H^, lowest in the case of T-D^ and intermediate 
in the case of D-H^. It is expected that the polarization potential arising out of a ls+2s+2p 
calculation will have the largest effect on convergence when the final-state velocity is the 
lowest. Hence the necessity of the higher-order states is more pronounced in the case of 
T-D^ and less pronounced in the case of T-H^. We also find that as energy decreases the 
transfer cross sections increase and the transfer rates attain a constant value. These transfer 
rates are essentially constant below 0.1 eV and are also measured experimentally, so that 
we can compare our rates with other experimental and theoretical results. 

For the D-H^ system the present low-energy muon transfer rate of 133x10^ s~^ is in 
agreement with both experiments [p!3| , p!4[ . The present rate is slightly smaller than the 



theoretical studies of Refs. pOl, |21] and this makes the agreement with experiment better. 



For the T-H^ system again the present result 61x10^ s ^ is in better agreement with the 



experiment |]T6| than the other theoretical studies. In case of T-D^, the present result 
2.3x10^ s~^ is also in very good agreement with experiment. 

Within the six-state approximation our cross sections for low energy elastic scattering 

10 



in case T-D^ system are presented in Table V together with other theoretical results. The 
present cross sections attain a constant value at low energies and is in fairly good agreement 
with results of other studies. 

As a futher test of the present few-body approach, we have also calculated S-wave cross 
sections of antihydrogen formation in antiproton-positronium low energy collisions (^. In 
Table VI our results within the six state approximation (Ps[ls+2s+2p],II[ls+2s+2p]) are 
compared with calculations based on hyperspherical coupled-channel method |26|. Consid- 
ering that the present calculation is limited to only the lowest partial wave {L = 0) and 
to a truncated basis set (ls+2s+2p), the agreement is reasonable for energies below 1 eV. 
However, at 2 eV the agreement is not so good. The reason for this is not clear at present. 
Further theoretical investigation including higher partial waves with an extended basis set 
could reveal the trend of the converged cross sections. 

IV. CONCLUSION 

The study of three-body Coulombic systems have been the subject of this work. We have 
formulated a method for a few-body description of the rearrangement scattering problem by 
solving the Faddeev-Hahn-type equations in coordinate space. It is shown that within this 
formalism, the application of a close-coupling-type ansatz leads to satisfactory results al- 
ready in low-order approximations for (i) muon-transfer reactions between hydrogen isotopes 
and (ii) antihydrogen formation in antiproton-positronium collision. Because of computa- 
tion difficulties, in this preliminary application we have considered up to six states in the 
expansion scheme (ls-|-2s-|-2p on each center), which may not always be adequate. Further 
calculations with larger basis sets are needed to obtain the converged results. 

The present model leads to a reduction of the usual technical effort and is definitely 
worth using for investigations of larger systems. It seems reasonable to suppose that the 
method should be an effective tool for the description of other muonic and atomic few- 
body collisions. For instance, one could study using the present approach the following 
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muon-transfer reactions to elements with Z >2 

(H^)i, + X^ ^ Xj + H, (29) 

where the cross section depends in a comphcated manner on the charge Z 0. 

Theoretically, the reaction ( PP| ) is of much interest as an example of low-energy rear- 
rangement scattering in a system of three charged particles with Coulomb repulsion in the 
final state. Evidently it makes additional difficulties for correct theoretical description of 
Eq. (p9D [|T^j . The Faddeev-Hahn-type approach seems to be suitable for the study of such 
reactions and would be a topic of future investigation. We are presently in the process of 
studying reaction ( ]29|) with the present method for Z = 2 and 3. We also plan to employ 
an extended basis set with more basis functions in the future. Also, the excited state muon- 
transfer reactions of recent experimantal and theoretical interest p9| , |30[| could be studied 
with the present model. 
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TABLES 
TABLE L Experimental and theoretical results for the muonic transfer rates A|;]-/10®s~^ given 
for low energies {E < 0.1 eV); *— present results; "— rates reproduced from cross sections. 



Reaction 



Experiment 



Theory 



D + H^ -^Bf, + H 



95 ± 34 H 

143± 13 

84± 13 ll| 



140 [go I 
133 [*1 



159 m 



T + H^^T^ + H 



58.6± 10 m 



55 [gO| 



61 



71.7 m 



T + D„^T„ + D 



M ^M 



2.9± 0.4 131 
2.8± 0.3 m 
2.8± 0.5 m 
3.5± 0.5 m 



3.5 [|0| 
1.5" (Hi 

2.39" m 



2.3 



2.26 m 
2.8 m 

0.93" H 



TABLE IL Cross sections f7(D-H^) = o-tr/lO^^'^ cm^ and rates A(D-H^) = Atr/lO^'^s^^ for 
^-transfer reaction D + H„ ^D„ + H, at different energies. 



E (eV) 


ct(D-H^) 


A(D-H^) 


a(D-H^) 


A(D-H^) 


^(D-H^) 


A(D-H^) 






Is 


ls+2s 


ls+2s+2p 


0.001 


292.6 


0.64 


412.8 


0.91 


604.8 


1.33 


0.01 


92.3 


0.64 


130.0 


0.90 


190.0 


1.32 


0.04 


46.0 


0.64 


64.7 


0.90 


94.3 


1.31 


0.1 


29.0 


0.64 


40.8 


0.90 


59.4 


1.31 


1.0 


9.0 


0.63 


12.8 


0.90 


19.4 


1.30 
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TABLE III. Cross sections a{T-B^) = a^^/W-^^ cm^ and rates A(T-H^) = Atr/lO^^s"^ for 
^-transfer reaction T + H„ -^ T„ + H, at different energies. 



E{eY) 


^(T-H^) 


A(T-H^) 


^(T-H^) 


A(T-H^) 


^(T-H^) 


A(T-H^) 






Is 


ls+2s 


ls+2s+2p 


0.001 


204.2 


0.42 


249.4 


0.52 


294.4 


0.61 


0.01 


64.3 


0.42 


78.5 


0.51 


92.6 


0.60 


0.04 


31.9 


0.42 


38.9 


0.51 


45.8 


0.60 


0.1 


19.9 


0.41 


24.3 


0.50 


28.6 


0.60 


1.0 


5.50 


0.36 


6.70 


0.44 


8.0 


0.52 



TABLE IV. Cross sections cr(T-D^) = atr/lO^^^ cm^ and rates A(T-D^) = Atr/lO^s"! for 
/i-transfer reaction T + D^ -^ T^ + D, at different energies. 



E{eY) 


^(T-D^) 


A(T-D^) 


a(T-D, 


) 


A(T-D^) 


^(T-D^) 


A(T-D^) 




1 


s 




ls+2s 


ls+2s+2p 


0.001 


4.58 


0.77 


5.05 




0.84 


13.7 


2.3 


0.01 


1.44 


0.76 


1.60 




0.84 


4.3 


2.3 


0.04 


0.71 


0.75 


0.78 




0.83 


2.14 


2.26 


0.1 


0.44 


0.73 


0.48 




0.81 


1.32 


2.21 


1.0 


0.1 


0.44 


0.1 




0.5 


0.3 


1.5 
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TABLE V. Elastic cross sections for T-D„ collision in units of 10 ^° cm^ at different energies. 



E{eV) 



Present results 
ls+2s+2p 



m 



0.001 
0.01 



1.2 
1.3 



1.7 
2.3 



1.63 
2.15 



2.014 
3.605 



TABLE VI. Cross sections in unit ttoq for positron transfer reaction p + Ps — > H + e" 
the cross sections estimated from Figure 1 of Ref . [26] . 



-• t_ 



E{eV) 

0.1 
0.5 



Present results 

ls+2s+2p 

1.5 

2.0 



26| 



2.3t 
2.6^ 



E{eV) 

1.0 
2.0 



Present results 

ls+2s+2p 

3.2 

1.7 



m 



3.5t 

3.7t 



